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Abstract 

We consider an optical diffraction grating in which the spatial distribution 
of open slits forms a fractal set. The Fraunhofer diffraction patterns through 
the fractal grating are obtained analytically for the simplest triad Cantor 
type and its generalized version. The resulting interference patterns exhibit 
characteristics of the original fractals and their scaling properties. 

(47.53. +n, 04.30.Nk, 01.55.+b) 

I. INTRODUCTION 

A widespread interest in fractal geometry has been generated by Mandelbrot with his 
monumental work |T| that deals with the geometry of phenomena observed in many fields of 
science. The fractal geometry is a fast growing subject and a cross-disciplinary field ranging 
from social science to biological science and physics |@-[|]- 

In this article we consider another aspect of fractal geometry, namely, the Fraunhofer 
diffraction of an optical grating whose slit distribution is fractal. In the ordinary optical 
gratings the slits are evenly distributed forming a one-dimensional lattice ||. We take, 
instead, a fractal grating and study the optical interference pattern by the light wave coming 
from the grating. As a simple example a Cantor set type grating is considered, and analytic 
results are obtained. The resulting interference patterns show self-similar structure and the 
same scaling property as those in the fractal grating itself. In other words the Fraunhofer 
diffraction transforms the fractal geometry of the optical grating to the fractal interference 
pattern. 

In Sec. II the simplest triad Cantor set type grating is introduced, and the interference 
amplitude is explicitly evaluated. The amplitude turns out to be a repeated product of a 
function of a rescaled argument by 3 m for each integer m, which respects the scaling property 
of the diffraction grating. In Sec. Ill we generalize the triad Cantor set to a (2H + l)-piece 
Cantor set, and similar results are obtained but with a different scale factor 2H + 1. In Sec. 
IV we discuss further generalizations of the model. The ordinary diffraction grating has been 
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an important tool in optical spectroscopy. Although we have not yet proposed any practical 
application, this kind of fractal diffraction grating may be useful for some purposes. 



II. CANTOR DIFFRACTION GRATING 

One of the simplest example of fractal sets is the Cantor set defined on a closed interval 
[0, 1]. The first stage of construction consists of dividing the interval [0, 1] into three equal 
pieces, then removing the middle open interval, designated (|, |). At the second stage one 
divides each remaining pieces into three and again remove each middle piece. Repeating this 
process infinitely many times one obtains a fractal set whose fractal dimension is D = 



Let us consider a diffraction grating whose slits are distributed as the Cantor set described 
above. For convenience, we take the starting total transparent region to be [—a, a]. At the 
first stage we make the middle third (— | , | ) opaque, in the second stage we again make the 
middle pieces of each of the two remaining thirds opaque, and so on to N times. The resulting 
grating consists of fractally distributed 2 Ar -slits of width J^-. Let this be called the Cantor 
diffraction grating of A-th generation. Figure 1 illustrates such a construction for the case 
of N = 3. The open slits are of equal widths but unevenly distributed. The coordinates of 
the middle, beginning, and end points of the j-th slit are denoted by Xj, xj , x~j~, respectively. 

The Fraunhofer interference pattern of light through a Cantor diffraction grating of A-th 
generation can be computed as 



/oo 
dxe lkxs[ne G(x) 
-oo 



2 ikx sin 9 + 

= V- \ Xj (l) 

fr[ ik sin 6* x s ' V ; 

where 9 is the angle of light propagation with respect to the normal vector of the grating 
plane, and k is the wave vector. The grating function G(x) is defined as 

r>( \ — / 0' x e opaque region, , , 

1 1, x £ open region. 

Using the fact that the slit width (x+ — xj) is J$-, we can write An as 

2 N 



2sin(^p^) 



A N = , v \ 1 Y e ikx > sine . (3) 

sin /feasm9j 

The amplitude An is factorized into two factors. The first factor, 2 — ^, si ^ e ' , is just the 
Fraunhofer diffraction amplitude of the single slit of width The second factor represents 
the interference of the 2 Ar -Cantor slits, which are expected to show the fractal nature of the 
slit distribution. 

In order to evaluate the summation in the second factor we first notice that the slit- 
coordinates Xj (the midpoint of the j-th open interval) can be represented by tertiary decimal 
numbers where only ±1 among (—1,0, +1) appears as numerator. Explicitly, they are 
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f-l -1 -1 -1 

Xl = 2a(— + — + -.. + — + — 

/-l -1 -1 +1 

x 2 = 2a[— + _ + ... + _ + — 

f -l -1 +1 -1 

x 3 = 2a[— + — + ••• + — + — 



^ = 2a( T + — + •■■ + ^ + —) . (4) 



Using these coordinates we can evaluate the summation as 



2^ TV 



^gttoyrinfl = -Q 2cO S ( 2fca f ng ). (5) 
j=l m=l 3"' 

Substituting to (Ef) we obtain the amplitude of the Cantor diffraction grating as 



2 ,N n sin(^) n /2 £;asin# 



"3 3 iV m= i 

= (1)^20^09)^09), (6) 
(i\ N 

where (3 = 2kasm9. Here (g) 2a signifies that the peak amplitude diminishes by the 

factor (|) as the open length of the grating decreases by the same factor. The Sn(/3) 
is the Fraunhofer diffraction amplitude of the single slit of the width J$-, and the F^((3) 
codifies the structure of the iV-th generation cantor set. 

The repeated product of cosine functions in the Fn(8) has a recursive nature as 

F^ +1 03)=cos(|)^(|), (7) 



In the limit N — > oo we have 



Foo(/9) = cos(^)F 00 (^) ) (8) 



which reveals the self-similarity with respect to the scale transformation by the factor 3. 
The maximum of F(j3) occurs at (3 = 0, F^(0) = 1, and the function oscillates very rapidly. 
Figure (2a) shows F N (f3) in the case 2ka = 3 10 , N = 10 for < 9 < §, while Figure (2b) 
shows the same function only for the region < 9 < 0.04. Apparently they look similar, 
and this self-similarity becomes more striking as N increases. The distribution of zeroes 
looks random, which may imply that the distribution of the dark fringes in the interference 
pattern is fractal. 



III. GENERALIZATION OF CANTOR DIFFRACTION GRATING 

For the purpose of generalizing the Cantor grating we first consider a function represen- 
tation of the Cantor set. Let g(x) be a periodic function of period 2a, 
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1, -a <x<-§, 
g{x) = { 0, -f<x<|, (9) 
-1, f < x < a, 



and define a grating structure function Gn(x) as 



N-1 



G N (x) = g(x)g(3x) ■ --g^x) = J] 9(^), (~a < x < a). (10) 



m=0 



Using this grating structure function the interference pattern through the Cantor grating 
can be represented as 



A N {(3) = P dxe ikxsin6 G N (x) 

J —a 



(l) N 2aS N (9)F N ((3), (11) 



where (3 = 2ka sin 9. In this way we construct a transformation between the grating function 
Gn(x) and the interference pattern function F N ([3) modulo the single-slit factor Sn(0): 

N-1 N a N-1 omn 

Gn{x) = J] 9(3 m x) ~ F N {(3) = n cos(^) = n cos(-f ), (12) 

m=0 m=l m=0 

which may be called a diffraction transformation between G^(x) and F^{f3). Their similarity 
in the product form of periodic functions and scaling property are obvious. To obtain the 
interference pattern we simply replace the generating function g{x) by the periodic function 
cos(^r). 

Let us now consider the possible generalization of the fractal structure. One immediate 
generalization is to take the generating function g(x) different from the simplest Cantor case. 
For example, we consider a periodic function h(x) of period 2a. As shown in the Figure 
3, the region [—a, a] is divided into 2H + 1 pieces of equal length with H being a positive 
integer. Its value is zero in the dark region and one in the open region, where the dark and 
open regions alternate. 

fl -g _|_ 4ma < < _ i (4m+2)a 
U(„\ — J ' ^ 2H+1 — X — U ^ 2H+1 ' (-\Q\ 

() [0, -«+<Sf <x<-a + ^, (m = 0,±l,±2,...). {U) 

With H — 1, one finds that h(x) reduces to g(x) of the previous section. From this periodic 
function we construct a grating function 

N-1 

G N (x) = J] h((2H + l) m x), -a<x<a, (14) 

which, in the limit — > oo, becomes a fractal set with the fractal dimension, 

\n(H+l) 

D= ln(2H + iy (15) 
The Fraunhofer diffraction amplitude through this grating is 
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O <j\ n ( kasme \ (H+1) N 

a n = y 2H+ » N) v e**>*«, (i6) 

A; sin [ 



where Xj are the midpoints of the open slits. Defining 



(H+1) N 

F N = e ikx i sind , (17) 



we find the recursive relation 

fn(p) = , Y 2H+ ; )N { f n -m, (is) 

S111 \(2H+1) N ) 

where (3 = 2ka sin 6. Proceeding as in the previous section, we finally get 

n-i sinf (H+1)/3 "l 

KM = n ) (2 T v ( 19 ) 

m =o sm( (2//+1)m j 
In the simplest case H = 1, we confirm the previous result 

n sin(|0 o 

*w) = n-fH = nscos a . (20) 

m =ism(^-J m= i x o 

The diffraction transformation then becomes 

n-i n s i n (AM+M.) 

Gn(x) = II h((2H + 1)™*) <- FtfGf?) = II , g±g^ . (21) 

m=0 m=l Sln ((2.f/-fl)™ ) 

We note that the factor function in F N is the familiar Fraunhofer diffraction pattern of 
a regular grating with (H + l)-slits whose width are / 2 H+i) m ■ This dearly exhibits the 
self-similarity with respect to the scale transformation by a factor 2H + 1. 



IV. DISCUSSION 

We have considered the Fraunhofer interference pattern obtained from diffraction grat- 
ings whose slits are fractally distributed. The simplest Cantor set type grating and its 
generalization were explicitly constructed, and their interference patterns were studied. The 
grating structure function, Gjy(x) = rim=o5 , ((2^+ l) m x^j, gives the interference pattern 

function, Fn(/3) = IIm=o /((zff+TV")' by the diffraction transformations, both of which show 
the repeated product structure, and have the self similarity with respect to the scaling by 
a factor of 2H + 1. However, the fractal dimension of the function of interference pattern 
F N (f3) is not fully understood. 

For general gratings with fractal distribution, analytic results may not be easily obtained. 
For example, if one tries the triadic Cantor generator whose segments are of unequal lengths, 
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the computations are no longer straightforward because the factorization is not possible. In 
this case numerical calculations would be interesting. If we consider a generator function 
g(x) of general form, it is not easy to picturize the product function G^(x) nor to evaluate 
the diffraction pattern function F N (f3). 

Another direction of generalization is to consider higher dimensional fractal diffraction 
gratings which we intend to study in future. 
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V. FIGURE CAPTION 

FIG. 1. Cantor Diffraction Grating for N = 3. The width of each open slit is ||. 
Xj,x~,x~j~ are the middle, beginning, and end points, respectively, of the j-th open slit. 

FIG. 2a. The amplitude factor F N (x) with x = 2kasm9 for N = 10, 2ka = 3 10 , and 
< 9 < f. 

FIG. 2b. The same amplitude factor F N (x), rescaled to the range < 9 < 0.04. The self 
similarity of the function would be more striking if N becomes larger. 

FIG. 3. The Cantor-like generator of (2H + l)-slit. It is a periodic function of period 
2a. The region [—a, a] is divided into 2H + 1 pieces of equal length. The h(x) takes values 
either zero or one alternatively. 
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